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We study the transfer of color between a heavy quark pair and an unpaired heavy quark or 
antiquark moving at a nonrelativistic velocity with respect to the pair. We find that the open 
heavy quark or antiquark can catalyze the transformation of the pair from octet representation at 
short distances to singlet at long distances. This process is infrared sensitive in general, and we 
exhibit double poles in dimensional regularization at next-to-next-to-leading order in the transition 
probability. Because of their dependence on kinematic variables, these poles cannot be matched to 
the non-perturbative matrix elements of effective field theories based on a single heavy quark pair. 

PACS numbers: 12.38.Bx, 12. 39. St, 13.87.Fh, 14.40Gx 

I. INTRODUCTION 

Heavy quarkonium production is of special interest in quantum chromodynamics (QCD), because the formation 
of the heavy quark pair leading up to its hadronization can be perturbative. The subsequent evolution of the pair 
has been extensively and fruitfully treated in the language of effective theories. Prominent among these are non- 
relativistic QCD (NRQCD) [H, 0, S, i and its extensions @, H, 0. In NRQCD, the non-perturbative dynamics is 
organized through matrix elements of operators that are characterized by an expansion in v, the relative velocity of 
the pair. 

The NRQCD expansion in v is well- understood for heavy-quarkonium decays although it still lacks a fully 
compelling proof for production processes. The NRQCD expansion, however, should not be expected to apply to 
reactions in which three heavy particles, two quarks and an antiquark for example, are produced close together in 
phase space, because in such cases there is no unique choice of relative velocity v. In this paper, we will extend the 
reasoning of Ref. Q, and argue that in this limited region of phase space, the formation of color-singlet quark pairs 
may be enhanced by a process that we will refer to below as color transfer. 

The application of NRQCD to production processes has had many successes [H, H, 0, H- There are, however, 
problematic cases. One involves the polarization of high-p^ charmonia as observed at the Tevatron [lol [Tlj . In 
connection with these data, the importance of associated production within the formalism of leading order (LO) 
NRQCD has recently been shown in Ref. [12J. Another unexpected observation was the large production cross section 
for J I tp in association with cc pairs, as seen by the BELLE and BAB AR collaborations [13[ . It has been shown that 
much of this excess may be accounted for within NRQCD by computing next-to-leading order (NLO) corrections [l4l |. 

The discussion below will explore the possibility that a complete picture for associated production will require us 
to extend NRQCD itself. We will argue that associated production with J/tp and related quarkonium states may be 
sensitive to the kinematic region mentioned above, where heavy quarkonium is produced close in phase space to open 
heavy flavor. The additional quark, itself a source of color, could influence the hadronization of these states. This 
mechanism also sheds light on NRQCD factorization [l5[ for processes without associated heavy quark pairs, and we 
will make contact with our previous work on that important question [l6l . |j"7| . 

In the following section, we begin with a discussion of infrared poles in dimensional regularization for NLO correc- 
tions to associated production in e + e _ annihilation We will exhibit infrared divergences that, although purely 
imaginary at this order, do not match with NRQCD operator matrix elements. These infrared poles are associated with 
the transfer of color between two heavy quarks and an antiquark, and cannot be associated with a quarkonium wave 
function. In Sec. 3, we analyze these NLO corrections in the limit that the relative velocity for one quark-antiquark 
combination is much smaller than that of the other, while both velocities are non-relativistic. In this non-relativistic, 
velocity-ordered region, we will exhibit an infrared divergence with a characteristic kinematic enhancement. We will 
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go on to identify the effective nonlocal operator that generates the leading contribution, and relate it to a similar 
combination encountered in the effective theory known as potential NRQCD [5j. 

Section 4 treats this color-transfer effect at next-to-next-to-leading order (NNLO) in the velocity-ordered region. 
We find in this case a real contribution to the cross section with a double infrared pole. Generalizing this purely 
perturbative result, in Sec. 5 we motivate an order-of-magnitude estimate for the color-transfer contribution to the cross 
section. We provide numerical estimates at B factory energies, which suggest that color transfer can be competitive 
with the familiar color singlet mechanism for associated production. 

We go on in Section 6 to study briefly the high energy behavior of the color transfer mechanism. Although we 
use the example of e + e _ annihilation, many of our observations apply as well to fragmentation in hadronic collisions. 
We exhibit why, as already found in [l2j], associated production is leading power in transverse momentum for the 
low-relative velocity region that we have in mind. We also touch on issues of polarization at high energy. We conclude 
with a summary that suggests directions for future work. 



II. NRQCD MATCHING AND NLO CORRECTIONS 



Nonrelativistic QCD organizes production and decay probabilities in terms of an expansion in the relative velocity, 
v <C f , of the quark-antiquark pair that eventually forms the bound state. In associated production, final states have 
two pairs. We will refer to the pair that eventually forms the bound state as the "active" pair. We will refer to the 
additional heavy particle that is closest in phase space to the active pair as the "spectator" quark (or antiquark). 

The quark and antiquark momenta of the active pair are Pi and P 2 , respectively. These momenta are related to 
the total and relative momenta and their rest-frame relative velocity by 

P P „ q 2 Am 2 

Pi = — + Q, P 2 = q, v 2 = -^ = \ (1) 

2 2 E* 2 P 2 y ' 

Here 2E* is the total energy of the active pair in the pair rest frame, and the second expression for v gives its relation 
to the quark mass and the invariant mass of the pair. 

Our discussion in this section will lead to a direct evaluation of familiar one-loop diagrams in the eikonal approx- 
imation, which will be followed in a later section by an analysis of the infrared behavior of all relevant two-loop 
corrections. In essence, we will supplement below the calculations of Refs. 0, \Yq . by studying infrared poles in 
corrections to the production of a heavy quark pair due to the exchange of soft gluons with an additional massive 
final-state particle. The extension from massless to massive particles in the final state has important consequences 
for a cross section in which the color of the quark pair is fixed to be in a singlet configuration in the final state. 
When other massive particles are observed, double infrared poles in dimensional regularization that are absent in 
the massless case appear in NRQCD coefficient functions. These poles have residues that depend sensitively on the 
relative momenta of the pair and the heavy particles. They are sharply peaked at small relative momentum, but fall 
off very rapidly when the relative momentum exceeds the masses. 



A. Lowest order and NRQCD 



Figure Q] shows the lowest-order diagrams that contribute to associated production in e + e~ — > QQ+QQ annihilation 
through a virtual photon. To form a heavy quarkonium state, it is assumed that the members of the active pair are 
close together in phase space, with relative velocity v <C 1 in Eq. (JTJ). 

To derive the NRQCD expression for the production of heavy quarkonium H corresponding to the diagrams of 
Fig. [TJ we first project these amplitudes onto a complete set of spin, angular momentum and color states, collectively 
labeled n, for the pair. We then square these amplitudes, normalize appropriately, and multiply each of the resulting 
expressions by the corresponding NRQCD matrix elements, which summarize nonperturbative dynamics. Summing 
over ro, we have, schematically, 

da e+e -^ H+x ( PH ) = da ^-^ Q Q[n ]+ x(PH) <0| Of |0) , (2) 

n 

with ph the quarkonium momentum. Color singlet NRQCD operators are of the general form [l| 

°5i)(°) = E X f (0)«n(i)^(0) \N,H) (N,H\ ^(0)< (1)X (0) 
JV 

= ^(QKajiKo) ^(o)< (1 )x(o), (3) 
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with K n and n' n projections for color (octet or singlet) and spin. For this factorization to be useful, it is necessary 
that the "coefficient function" da e + e -^QQ^ +x (pH) be infrared safe at higher orders. This factorization remains a 
conjecture for production processes beyond NNLO. 

In Ref. [HI] we showed that, beyond NLO, factorization requires a reformulation of octet production matrix elements 
into a "gauge-completed" form, 



n(8)(0) = x t (0) K „(8)x^(0)<i>|' 4)t (oo,0) cfc (a^an) 3>\ A) (<x>, 0) ba ^(0)< (8) ,j(0) : 



O 



$[ A) (Ai,0) = V exp 



-ig d\l^A^ a (Xl)T^ 
Jo 



(4) 



Here, as shown, ^ A \\, 0) a h is a gauge link in the adjoint representation (a, b = 1 ... 8) with a light-like four- velocity 
I, and with V path ordering for the expansion in g. We showed, in particular, that gauge-completed matrix elements 
are independent of the direction of the gauge link, so long as I remains light-like. We will see in the calculations that 
follow, however, that independence of the direction of / does not generalize to a massive color source on the right-hand 
side of Eq. ([JJ . This observation will have important consequences for associated production. To see how this comes 
about, we turn now to soft gluon loop corrections to the lowest order diagrams of Fig. [T] 



B. NLO soft gluon corrections and matching 



It would be a major undertaking to compute the full NLO corrections to the singlet NRQCD cross section for two 
pairs at arbitrary momenta. Nevertheless, it is relatively straightforward to check the self-consistency of the NRQCD 
factorization at NLO in this process. This can be done by checking that infrared poles in dimensional regularization 
either cancel, or can be absorbed into vacuum expectation values (O|0 n |O), thus matching full QCD to NRQCD. 
Matching is essentially equivalent to NRQCD factorization. 

More specifically, NRQCD factorization and matching for production cross sections require the cancellation of 
all infrared gluons that are not "topologically factorized" into factors equivalent to the perturbative expansions of 
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matrix elements in the effective theory [lj, The allowed, topologically factorized, soft gluons are those that can 
be absorbed either into the interaction of the active pair with "the vacuum" , or that do not couple to either of the 
heavy quarks that form the quarkonium. In the former case, the soft gluons generically will have the interpretation of 
part of a nonperturbative matrix element or wave function. In the latter case, they will cancel in the inclusive cross 
section for fixed active pair color representation. 

Figure [2] illustrates these considerations. With Pi and P2 the momenta of the active pair, the diagram on the upper 
left of the figure shows a virtual soft gluon that is emitted from the active quark and absorbed by the active antiquark. 
Such a correction can be absorbed into a wave function. Similarly, contributions that describe the interference between 
gluon emission by the active quark in the amplitude with emission by the active antiquark in the complex conjugate 
amplitude are also topologically factorized, and hence consistent with NRQCD matching. 

The right of the upper line of Fig. [2] shows a gluon that connects to neither of the active quarks. Such a contribution 
cancels in the sum over final states at fixed color representation of the active quarks. It is thus also consistent with 
NRQCD matching. 

In the remaining two diagrams of Fig. [2j a member of the pair interacts with a spectator quark, of momentum P3 in 
this case. These two diagrams are therefore not topologically factorized. Nevertheless, in a production cross section 
the real infrared poles of these diagrams still cancel against diagrams with soft gluons in the final state, for a fixed 
color projection on the active quark pair 'lj. Of course any imaginary terms cancel upon combination with complex 
conjugate diagrams. At NLO, therefore, the non-topologically factored diagrams do not affect NRQCD calculations 

In Ref. [17J, however, we found a somewhat stronger result for the imaginary poles of the non-factored diagrams of 
Fig. [21 When the spectator momentum is light-like (i.e., P| = 0), the "Coulomb phase" associated with the exchange 
of a soft gluon between a massive quark and a massless quark is independent of their relative directions. As a result, 
the imaginary infrared poles of the two non-factored diagrams in Fig. [5] differ only by a relative minus sign between 
the quark and antiquark in the active pair. Once we project on a singlet final color state for the pair, they then cancel 
identically, to all orders in the relative velocity, v, Eq. ([1]). This is an important ingredient in the NNLO factorization 
of gauge-completed octet NRQCD matrix elements, Eq. (QJ. 

We will reproduce these diagrammatic results below, but show that they no longer hold for a massive spectator 
(anti)quark. For a massive spectator, the infrared pole is a simple but nontrivial function of the pair relative velocity 
v, as well as the pair-spectator relative velocity. 



+ 




FIG. 2: Diagrams for the one-loop virtual infrared corrections, where the blob represents a complete set of lowest order Feynman 
diagrams, as shown in Fig. [1] 

At the level of infrared poles in dimensional regularization, the one-loop virtual corrections for the diagrams shown 
in Fig. [5] can be factorized into LO short-distance color matrices, M( LO ), times long-distance color matrices, which 
describe the exchange of soft gluons. The latter can be calculated in the eikonal approximation, in which the jth 
quark or antiquark propagator and vertex are replaced by the spin-indendent combination 




(5) 
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with the plus sign for a quark and the minus for an antiquark, and where / = q or q denotes the quark or anti- 
quark representation for the generators Equivalently, we may take all color matrices in the defining (quark) 
representation, and stipulate that matrix multiplication follows the arrows that represent the flow of fermion number. 
In these terms we write the factorized amplitude as 

4T' IR) = E ^^^Uii^^^^Ps^^M^CPr^s^a^^, (6) 

where Af ( LO ) is the leading-order amplitude projected on the appropriate spin state, at fixed values of color indices. 
A^ 1 ' is the eikonal factor, describing the coupling of a soft gluon to the outgoing quarks. The diagrams of Fig. [3] are 
typical contributions to A^\ and in this case correspond to the non-factored diagrams in Fig. [2J with the heavy lines 
representing the eikonal approximation ([5]) for the couplings of soft gluons of momentum k. 

In Eq. ©, color indices i\ and ii correspond to the active pair, of momenta Pi and P2, which are associated with 
the heavy quarkonium. As shown, we set i\ — 12 and sum to enforce a color singlet final state. The remaining color 
indices 13 and 14 correspond to the open heavy quark pair, of momenta P3 and P4, as shown in Fig.[TJ We now turn 
to the evaluation of the non- factored contributions. 




FIG. 3: Diagrams with eikonal interactions between the heavy quark pair and an associated heavy quark of momentum P3. 
The heavy lines indicate fermions in the eikonal approximation. 



C. Infrared poles in the NLO amplitude 



The contributions of the two diagrams in Fig. [3] to the eikonal factor in Eq. © are given in D = 4— 2e dimensions 

by 



d D k p 1 ■ p 3 



(2tt) d [Pi • fc + ie] [-P 3 -k + ie] [k 2 + ie] ' 

d D k P 2 -P 3 

(2n) D [P 2 -k + ie] [-P 3 -k + ie) [k 2 + ie] ' 



(7) 



with A = A13 + A23 + . . . . For the analysis of this section, we will assume equal masses m — P 2 = yPf for the 
active quark pair, but allow the additional heavy flavor to have a possibly different but nonzero mass, 7713 = v/Pf . 



3 • 

In accordance with our previous discussion, we will be interested in the relative motion of the three heavy particles. 
For the equal-mass case, a convenient choice for relative velocity between Pi and Pj is /3y = -jT— Am 2 / Sij, where we 
adopt the notation 

% = (P + Pi) 2 . (8) 

As in Eq. ([1]), we define v = P12 for the relative velocity between the active pair. Consider now momenta p and p' 
with arbitrary masses, m and m! . Their relative velocity in any frame where their spatial components are collinear is 
given by 

M = |v-vL = ^v¥7FHT, (9) 
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where E = po, and similarly for E' in the same frame. For equal masses and in the center of mass frame, (3(Pi, Pj) 
reduces to twice the familar relative velocity, y/1 — Am 2 /sij, of Eq. (P) . 

Another measure of the distance in phase space between two vectors that arises in the soft gluon corrections at 
hand is 



P(p,p') = ^P(p,p') 
P-P 



1- 



• p') 



(10) 



We easily verify the limiting behaviors, 



lim 0(p,p') = P(p,p') , 

(mm') /{p-p , )—*l 

lim P(p,p') = - P{p,p') ■ 

(mm' ) / (p-p' )—>0 Z 



(ii) 



That is, for small relative velocities, (3 is nearly equal to (3, although it increases somewhat more slowly with center- 
of-mass momenta, reaching a limit of (3 ~ 1 for fully relativistic motion. We will be interested in regions of phase 
space for which relative velocities are nonrelativistic. In such regions, /3,/3 <C 1. For convenience, we will refer to [3 
as a velocity. 

So that we may unambiguously distinguish the spectator quark, we restrict ourselves to momenta for which v < 
(3 « 1. In this region of phase space, conventional NRQCD factorization need not apply directly, because of the 
presence of an additional small parameter, (3. Indeed, we will find an additional infrared-sensitive enhancement in 
the octet-to-singlet conversion for the active quark pair, due to the presence of a nearby spectator. If this mechanism 
is observable, it should manifest itself in a peak in associated open heavy flavor distributions in phase space near 
quarkonia. 

The evaluation of each integral in Eq. is reasonably straightforward, and we readily isolate their infrared poles. 
In fact, in dimensional regularization these integrals vanish through the cancellation of IR and UV poles. Nevertheless, 
the single infrared pole may be found (for example) by adding quadratic terms +k 2 /2 to the eikonal denominators. 
The result is 



Pole( IR ) 



-19 



d D k 



P ■ P< 



(2ir) D [P t -k + ie] [-Pj ■ k + ie] [k 2 + ie] 

l+p(Pi,Pj) 



1 a* 



1 



2e 2wp(P i ,P j ) 



hi 



- 2in 



(12) 



where (3 has been defined in (fTTJ]) . In the region that we have identified above, /3(Pi, Pj) <C 1, we observe that for the 
infrared pole, singular behavior in /3 appears in the imaginary, but not in the real part, 



Pole( IR ) 



d D k Pi- Pj 

(2n) D [P t -k + ie] [-Pj ■ k + ie] [k 2 + ie] 



(13) 



Our discussion below will concentrate on the effect of this pole on production cross sections, which, as we shall show, 
is seen first at NNLO. 

At this point we note that the real part of the full vertex correction has a (famous) 1/(3 singularity. This power 
singularity near threshold, however, is finite in four dimensions, and is associated with loop momenta at the scale 
of iiro 2 , the so-called "ultrasoft" momenta. Dynamics at this scale are regulated by bound state effects between the 
active quarks. In this paper, we shall not attempt an analysis of exchanges at the ultrasoft scale between the active 
pair and spectatators. 

Consider now the imaginary contributions of Eq. J7]) to A, traced over the colors of the active pair to enforce a 
singlet configuration in the final state, 



U*2 



Im L4i 3 + A 



2:i_ 



l\...H\l, ...l. 



£ ( 2 ) ^ la ^ i 2 i i ^ a ^i3i' 3 ^ 4 *4 



y/1- PfPZ/fr.Ptf 



^1-P 2 P 2 /(P 2 -P 3 ) 2 



(14) 
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This pole in the imaginary part vanishes identically for P% ■ P3 = P2 • P3 or for P| = 0. In the former case, the 
relative velocity of the spectator quark and the active quark equals the relative velocity of the spectator and the 
active antiquark. The interpretation of the cancellation is then that soft gluons emitted by the spectator cannot 
resolve the charges of a perfectly co-moving pair in a singlet color state. Analogously when the spectator is lightlikc, 
the relative velocity to both active lines is unity and the two terms cancel. Evidently, in this case a light-like spectator 
cannot resolve a color singlet even when the quark and antiquark are not co-moving. 

From Eq. (114|) , the sum of the exchanges between the active pair and the spectator depends directly on the relative 
velocities of the spectator with the active quark and antiquark, but only indirectly on the relative velocity v of the 
active pair itself. In general, for /3i3 ~ /?23, we must describe the kinematics of two quarks vying for the favor of 
a single antiquark to form a heavy bound state, with no obvious favorite. This three-body problem is potentially 
difficult to analyze. The soft gluon dynamics of two quarks and an antiquark might be generated from local operators 
like *ipp 1 \p 2 Wp j with time evolution generated by heavy quark effective theory for each of the quarks, but only if we 
neglect the recoil of the quarks. Nonrelativistic QCD organizes such corrections into nonperturbative matrix elements, 
but only to the extent that the extra quark may be neglected. 

We will not attempt a full solution to the three-particle problem, but will concentrate in the following on a limit 
more closely related to NRQCD. We have already noted that the sum of NLO phases in Eq. (114|) vanishes in the limit 
v = 0i2 ~> 0, at fixed, finite values of P13 = 023- We now turn to this region of nonrelativistic, ordered velocities, 
which we will refer to as the "velocity-ordered region" . 



III. NLO COLOR TRANSFER FOR ORDERED VELOCITIES 



We begin this section by expanding the residue of the infrared pole on the right-hand side of Eq. (|14[) to lowest 
order in the active pair relative velocity, v, in the velocity-ordered region. The residue is strongly peaked toward low 
relative velocity for the spectator pair. We will then go on to rederive this result from an effective non-local vertex 
involving the gluon field strength coupled to the dipole moment of the active pair [l6| . This reformulation will enable 
us to extend our analysis to NNLO in the following section, and also to exhibit the relation of the analysis here to 
our previous study of NRQCD factorization [l7| and to potential NRQCD 



A. Expansion in v 

To estimate the kinematic behavior of the poles in Eq. (|14|) in the velocity-ordered region, we expand at low relative 
velocity, v for the active pair. In Eq. (114|) . the only kinematic variables are Pi ■ P3 and P2 ■ P3, and we expand both 
around their values at v = 0, (-P/2) ■ P3, with P 2 = 4m 2 We will allow P3 to have an arbitrary nonzero mass, 777,3. 
Given that we have two invariants, we will have to introduce an additional relative velocity, and corresponding relative 
momentum. 

We begin this straightforward kinematic analysis with Eq. (fl]) for the momenta of the active pair, recalling that 
the mass-shell conditions P 2 = P| = m 2 require q ■ P = 0. In a frame where the total pair momentum is at rest, the 
relative momentum, has no time component. We can then write the total energy of the pair as a function of q, 
in terms of a momentum at rest in that frame, Pq = (2m, 0). In covariant form the dependence of P on q is given 
explicitly by 



q 



2 



where q 2 = — <f c 2 m . We now parameterize P3" in a similar way, expanding around P3 = (777,3, 0) in the P M (active pair) 
rest frame, in terms of a vector g^, which like q has vanishing time component in that frame, 




^=^3^ a/1- 95-Po-O. (16) 



Expressing P3 in these terms, and recalling that P\ — P/2 + q, we evaluate the relevant invariant P\ ■ P3 as a function 
of q and qs, finding 



P-P 3 - 777 7773 J(l--C) f 1 -^) +1-1S- (17) 
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The analogous result for P2 • P3 is found by simply changing the sign of q ■ qs- We are now ready to expand (fT"4f in q 
at fixed qs- 

Starting at q = 0, P\ = P2 = Pq/2, both invariants P, • P3, i = 1,2 are given by 



(Po/2) • P 3 = m J ml -<%= -^L= , (18) 
where in the second expression we recall the definition of the function (3, Eq. (p~0|) . and define 



Ps = 0(Pb/2,ft) = J-JK. (19) 

V m 3 *?S 

This is a measure of the spectator's velocity relative to the active pair, neglecting the latter's internal relative veloctiy, 
v. We note that this particular velocity is independent of to, so that, for example, /3(Po/2, P3) = (3(Po, P3), and 0s = 1 
for m 3 = 0. 

Following our previous discussion, we consider the region q 2 < q s , and expand in the corresponding ratios on the 
right-hand side of Eq. HH), keeping in mind that at q = q$ = both terms in the difference on the right-hand side 
diverge. We note as well that only terms that are odd in q^ survive in the difference. We find 

1 1 2 IS ■ Q ( , 52\3/2 |\ , m ( q 2 



sjl - pJpJJWi ■ P3) 2 ^WpJKpFW P 3 s mm * 



1 + 



(20) 



This result is exact in qs up to terms that are quadratic in q. Notice that it vanishes rapidly in the limit m^/qs — > 0, 
that is, relativistic motion for the spectator. 

We can now define the velocity-ordered region as 

Ps<l, W <1 - ( 21 ) 

PS 

That is, we shall assume that it is possible to expand in v//3s, keeping in mind that not all v or fis dependence is of 
this form. Of course, such an expansion may be of limited quantitative use when j3s ~ v. This is a region likely to be 
of particular importance for charm quarks, where even v need not be small, and where the dynamics goes over into a 
true three-body problem. Nevertheless, we hope to gain insight even from the somewhat idealized kinematical limit 
of Eq. (ED). _ 

For small /3s, Eq. (fT9|) shows that the vector qs is proportional to /3s, just as the active pair's relative momentum, 
q is proportional to v. In this way, we can rewrite our result in a form that exhibits the leading dependence on both 
velocities v and (3s in the velocity-ordered region, 



1 1 2 A 

= — V COS <ps 



y/l-PfPf/fr-Ptf y/1- P*P*/(P 2 -P 3 )* (3 S 



771 771 TOq 



(22) 



where on the right 4>s is the angle between qs and q in the rest frame of the active pair. Equation (|22[) provides an 
explicit illustration of the enhancement of gluon exchange as the spectator approaches the active pair in phase space. 

We can summarize our analysis for the imaginary NLO pole in the amplitude at small qs by using (12TJ1) and (f2"2")l in 
(fT4)l . which gives 

» <•' 1111 A > ' - ""(if) ( Ta ki'A T ^' 3 ^ F (l-^) 3/2 «COS^ S . 

ii%2 S 

(23) 

We will return to this expression in Sec. 5, when we discuss its possible phcnomcnological implications for quarkonium 
production in association with open heavy flavor. In the remainder of this section, we rederive the basic result (|23|) 
by expanding first at the diagrammatic level. 
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B. Operator interpretation and potential NRQCD 



We have observed that the eikonal factor, A of (|T4[) vanishes at v = 0, that is, at vanishing active pair relative 
momentum. The lowest order in v exhibits an infrared singularity characteristic of an electric dipole transition [l6l |. 
although in this case the dipole is coupled to the field of a spectator rather than to an emitted gluon. The expansion 
of Eq. (|14p shows this result at lowest order, through its linear behavior on the relative momentum q. In this section, 
however, we will expand the diagrams first to order q, and sketch the resulting calculation. We may think of this as the 
first step in the creation of an effective theory for color transfer, in the restricted region where v < /3(Pi, P3), (3{P2,P?,). 






FIG. 4: Graphical representation of the first-order expansion in the relative velocity of the active pair. The vertex Cg> represents 
the inner product of the chromoelectric field and the relative velocity, as in Eq. 



We represent the sum of the two eikonal appoximations to the coupling of a soft gluon to the active pair in the NLO 
eikonal factor in Fig. 2J As in Ref. the double line on the right-hand side of the figure represents the combined 
quark pair. As we shall see, the double line and the vertex ® corresponds to a Wilson line in adjoint representation 
that terminates at the field strength operator P^F^ ' q vi where q is the relative momentum of Eq. |T]). 

Returning to the NLO integrals of Eq. 0, and recalling the relation between q^ and v in Eq. |T]), we see that an 
expansion in powers of v is equivalent to an expansion in q. For this purpose, we denote the spectator momentum by 
I, which may represent P3 or P4 in this case. This was the notation used in Refs. [H, where, however, we took 
l 2 = 0. 

In these terms, the order-g correction to the two contributions to the eikonal factor A of Eq. ^ is given by 

A 3 +^23 - M J {2n)D [ p. k + ie] 2 [ _ l . k + ie][k 2 + ie] + U \ P 2) ' W 

where the expansion results in a squared denominator for P ■ k and in the numerator factor 

N(P, q,l,k)=4[q-lP-k-P-lq-k] . (25) 

As shown in Ref. (lo| . we can interpret this expression as the local field-strength vertex appearing in the nonlocal 
operator 

/•oo 

E b = -ig dXX [?W(^)]^ ) (A,0) a6 . (26) 
Jo 

In the rest frame of the active pair, this is precisely the integral along a massive world-line, specified by momentum P, 
of the scalar product of the relative momentum with the electric field operator, contracted with a gauge link in adjoint 
representation (defined as in Eq. Q). Notice the factor A in the integral, which reflects the increasing separation of 
the quark and antiquark as the pair propagates into the final state at fixed relative velocity v = q/E* (Eq. |T])). The 
extra factor of A produces the squared denominator in Eq. |24|) . 

In Eq. (|26|) . the connection of our NLO correction to the effective theory potential NRQCD (pNRQCD) [5J] is 
manifest. The gauge link is the propagator for an octet pair in pNRQCD, and the field-strength is the operator in 
the Hamiltonian that transforms an octet pair to a singlet configuration. Although it was not designed to describe 
quarkonium production, pNRQCD operators emerge naturally in the velocity-ordered region, as a description of the 
active pair evolving in an external color field supplied by the spectator. 

We now turn to the evaluation of the NLO correction starting from the effective vertex of (|2"B)) . It is straightforward 
to evaluate the integral in Eq. ([24]) , but we describe it in some detail because it sheds light on the role of the spectator 
mass. For simplicity, we may choose l± = and P = (Mu, 0). From Eq. (|25| . we then find, after a little algebra, 

N = 4 M H q 3 [{I ■ k) - 21+ k-] +AP-lqj_-k ± , (27) 

where the last term vanishes after integration over k±. The first term in square brackets, I ■ k, cancels the spectator 
denominator. Because this term is actually independent of the spectator momentum, it can be treated in exactly the 
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same manner as in the massless case [17|, leading to the same result. This contribution is purely real, and cancels 
against real-gluon emission. We are left with the contribution of the l + k~ term only, which we now evaluate. 

As in Ref. [lB], we integrate over k~ first. When we close the k~ integration contour in the lower half-plane, we 
encounter poles from both the k 2 and P ■ k denominators. The k 2 pole gives a real result and cancels against the 
corresponding numerator factor for the final state with a gluon. The remaining contribution in the l + k~ term in Eq. 
([24]) is from the double pole at P ■ k = in Eq. ([24]) , which gives 
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k\ - ie] [V2l 3 k+ + ie] 



2{k+) 2 + k\ ^2P 



(28) 



In the second equality we have taken the derivative with respect to k~ . In the resulting expression on the right, the 
first term in brackets vanishes, because it is odd in k + , while the second term can be found from the pole at k + = 0, 
and has vanishing real part. After integration over k±, we obtain an imaginary piece for the eikonal factor, 



i = Im 
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= - 2a s 
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(29) 



found here in the P = (Mh,0) frame with £± = 0. In a general Lorentz frame, the leading order eikonal factor is 



Im 



A 1 



= -- (2a 



P 2 l 2 



[{P -l) 2 - P 2 l 2 f/< 



(l-Q) 



(30) 



This result, of course, agrees (up to color factors) with the expansion of Eq. (|14| for the imaginary parts of the 
diagrams in Fig. [3] As in Eq. ([H]) , the essential role of the spectator mass, VI 2 <-> is manifest. We also see that 

for P ■ I > I 2 , the effect decreases as (P ■ Z)~ 3/2 . 

Having arrived at a result identical to the expansion of Eq. (fT4|) at NLO, we are now ready to use the operator 
formalism to study color transfer at NNLO in the velocity-ordered region. 



IV. COLOR TRANSFER AT NNLO IN THE VELOCITY-ORDERED REGION. 



So far, we have identified an IR pole in the NLO imaginary part of the production amplitude for a color singlet 
(active) quark pair, starting from an octet pair at short distances, by using the eikonal approximation. This contri- 
bution depends on the presence of a massive source, the spectator, nearby in phase space. The spectator catalyzes 
the octet-to-singlet transition of the pair, in the sense that its own color representation is unchanged. The effect 
decreases rapidly as the relative velocity of the source approaches unity, where the amplitude is readily matched to 
NRQCD. On the other hand, when the relative velocity between the pair and spectator is itself nonrelativistic, the 
effect is potentially significant. 

Following the logic of Ref. [l6[ , we explore the implications of these infrared poles for the production rates of heavy 
quarkonia in association with open heavy flavor. We will work in the velocity-ordered region, (|21|) . where the active 
pair's relative velocity, v is even smaller than its velocity relative to the spectator. We then expand in v. 

From the point of view of NRQCD, we will calculate soft gluon corrections to the coefficient function 
d&A+B—>QQ[n]+x(PH) m ^1- ©' ^ or ^ ne specific case of QQ[n] a pair of heavy quarks in singlet configuration. We 
will identify real infrared double poles in this production amplitude starting at two loops. As above, at the level of 
infrared poles in dimensional regularization, these can be computed using eikonal approximation, or equivalently, as 
vacuum expectations values of products of Wilson lines, defined in adjoint and fundamental representations as in Eq. 

©■ 
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A. Matrix elements for pair production at NNLO 

To recall, Ref. [Tt| dealt specifically with infrared poles in fragmentation functions. The relevant eikonal production 
fragmentation function is given to lowest order in q ~ mv by 

X^ l) {P,q,e) = 2^ r d\'\' (0| $| A) t(oc, Q) bd , <^ A) (A', 0)^, [P^F v „, a ,{\'P)] \N) 

N J ° 

x(N\f dXX [FVF I , M , a (AP)]$^ ) (A,0) Qd $[ A) (oo,0) (ib |0), (31) 
Jo 

with (anti-) time ordering implicit in the (complex conjugate) amplitudes. Here $\ A ^ is the same light-like Wilson 
line in adjoint representation link as in Eq. (f4]), which represents the effect of recoiling massless quanta. The other 
Wilson line in adjoint representation, 5>p^, is massive, and represents the propagation of the pair of heavy quarks 
as an octet. This gauge link connects the hard scattering at the origin with the field strength tensor as in Eq. (I26[) . 
which describes the absorption of the soft gluon that changes the net color of the pair from octet to sing let. We note 
that in this paper we choose the argument of the field strength to be X'P, rather than X'P/2, as in [la ], which leads 
to an explicit factor of 2 on the right-hand side of (|3"Tj) . 

In this section, we study an eikonal production cross section, in which the light-like gauge link of Eq. (f3"Tj) is 
replaced by a pair of massive Wilson lines in quark and antiquark representation. These two lines are linked together 
by an octet short-distance vertex (the shaded circles of the figures), which matches to the adjoint representation of 
the line in the same direction as the total momentum, P, of the pair. Again, the adjoint Wilson line terminates at 
the field strength vertex, Eq. (|26|) . 

/•oo 

M 2 (P,q) ee ]T / dX'X' (0| $gt(oo,0)^ (T c ,) rk , 0>®t(oo,0) fc ^^ ) t(A',0) c w [P»q v F^ a ,(X'P)] \N) 

N Jo 

x (N\ / dX'X' [P»q»F^, a (X'P)]$¥\x\Q) ac $®(oc,0) lk {T c ) kj $g (00,0)^0) , (32) 
Jo 

again with (anti-) time ordering implicit in the (complex conjugate) amplitudes. This production cross section at 
NLO is represented by the diagram on the left in Fig. \5\ As in Fig. [4j the double line stands for the pair of heavy 
quarks at low relative velocity v, terminating in a singlet final state represented by the vertical line in the figure. The 
vertex (g) supplies the necessary color. We sum over all final states that include the singlet pair (P) and the antiquark 
line P4. 

The diagram on the left of Fig. Ogives the entire NLO contribution to A4 2 , Eq. (f32|) . It represents the classic 
lowest-order octet mechanism, in which the pair in octet representation emits a single gluon and becomes as singlet. 
It has a single infrared pole in dimensional regularization, and is the lowest-order contribution that is matched to 
nonperturbative matrix elements in NRQCD [19j . At this order, the associated pair is irrelevant to the infrared 
structure, so long as we concentrate on octet short-distance functions. 

The diagram of the same order on the right in Fig. [5] represents an interference between an octet pair produced at 
short distances in the amplitude with a singlet pair produced at short distances in the complex conjugate amplitude. 
In the amplitude, the pair is transformed into a singlet by the exchange of a soft gluon with the spectator quark 
line P3. We include this interference diagram because it illustrates at NLO the cancellation of infrared poles from 
an on-shell gluon (fc 2 = 0) in virtual corrections and in the final state. Both of these states are available even when 
the heavy active pair (double line in the figure) is fixed to be in a singlet state. After the cancellation of the real 
gluon contributions, however, the imaginary pole that we studied in section 3 above remains, and must be cancelled 
by complex conjugate diagrams (not shown in the figure). 

Although the pole we have found above at one loop is purely imaginary, it will contribute to the real part at NNLO. 
This pole, of course, is just the one-loop contribution to a Coulomb phase, and is guaranteed to cancel in a fully 
inclusive cross section. As we shall see, however, when restrictions are placed on the color of final-state pairs, as in 
NRQCD, this cancellation fails in general at NNLO. 

B. The analysis of double poles at NNLO 

The NNLO diagrams that contribute to the squared octet-to-singlet transition probability, Eq. (|32|) are shown 
in Fig. [SJ In these diagrams, the cut double line represents the active pair in color singlet configuration, and the 
crossed circles the field-strength vertices. All lines in these diagrams except for gluons are eikonal. The transition 
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FIG. 5: NLO contributions to 7* — > [QQ]QQ with the double line and the vertex g> defined in Fig. [4] The cut double lines in each 
case denote a heavy quark pair in a singlet state. The shaded circles of the diagrams represent short-distance functions in the 
amplitude and complex conjugate. In the figure on the left, both short-distance functions are octet, corresponding to the matrix 
element M2 of Eq. (I32p . On the right, the short-distance function of the amplitude is octet, and of the complex conjugate, 
singlet. The right-hand figure, which illustrates the cancellation of real and virtual on-shell gluons, does not correspond directly 
to Eq. (Eg) . 



probability (|32[) has been constructed to reproduce the infrared behavior of diagrams in full QCD at leading order in 
the active-pair relative velocity, v. 

The graphs shown in Fig.[B]each represent a set of cut diagrams, found by summing over all final states that include 
the color singlet pair and the line P4. As in Fig. [5j the shaded circles represent short-distance functions, which link 
the active pair in octet representation with the associated pair (P3 and P4). We work in the velocity-ordered region, 
Eq. (|2"Tj) . and take the line labeled P3 as the spectator that interacts with the pair. We will neglect connections of the 
active pair with the remaining line, P4. That is, we assume that the relative velocity of this line to the pair is large. 

Our goal in this section is to show that at NNLO there is a real infrared-sensitive contribution to the octet-to-singlet 
transition amplitude, M.2- We will not calculate the diagrams explicitly, but argue that in dimensional regularization 
diagram (I) produces a double pole that equals the absolute square of the NLO result, Eq. (|23[) (up to color factors), 
and that there is no other source of real, double infrared poles from any other diagram. Such a real contribution to 
the transition amplitude cannot be matched to a local operator in NRQCD. We will rely heavily on power counting 
analysis, and on the specific calculations of Ref. fl6| . 

To generate a double pole in D = 2 — e/2 dimensions, a diagram must contain two momentum integrals, either 
virtual or real, that diverge logarithmically when D = 4. This can be checked by straightforward power counting. 
General rules for power counting estimates of infrared-sensitive diagrams are discussed, for example, in Ref. [18j . The 
prescription is particularly simple for massive cikonal lines exchanging gluons without gluon self-coupling vertices. 
This describes all the diagrams of Fig. [5] except for (III). Because all the eikonal lines are massive, there are no collinear 
singularities anywhere in loop momentum or phase space for diagrams (I), (II) and (IV)- (VI). In contrast, we do find 
collinear singularities for some momentum configurations of diagram (III) because of the three-gluon coupling in that 
diagram. These collinear singularities, however, cancel in the sum over cuts at fixed active pair color by standard 
unitarity arguments [Tc| , and we need not consider them here. 

To perform power counting for the remaining, infrared singularities, we need to study all choices of loop momenta 
carried by the gluons. The momenta of the active pair is fixed, as are P3 and P4. When a gluon appears in the final 
state, we do power counting on its phase space integral, which we shall continue to refer to as a "loop" momentum. 
From a technical point of view, in the definition of M.2, the vertices associated with the short-distance functions 
are fixed at the origin in both amplitude and complex conjugate. These vertices are therefore relatively local, and a 
loop momentum can begin at the short-distance function in the amplitude on the left of a diagram, and end at the 
short-distance function on the right, in the complex conjugate amplitude. In physical terms, because we neglect the 
recoil of the heavy quarks, we do not impose overall momentum conservation on soft gluons in the final state. 

We scale the components of any such choice of loop momenta as 

K = Aifcf , (33) 

where i = 1, 2 labels a choice of loop momenta. This is just a change of variables, which will produce an overall factor 
A^ for each loop. When the integrand behaves as A~ Pl with pi = 4, all sizes of loop momenta contribute equally to 
the result, and the integral diverges logarithmically. In fact, because these eikonal integrals are scaleless overall, there 
is always a single logarithmic divergence in both infrared and ultraviolet when both loops are scaled together, and 
the diagrams are nonzero in dimensional regularization only after imposing an ultraviolet cutoff. 

The application of the scaling (|33[) is quite straightforward, keeping in mind several technical observations. 

1. To derive a double infrared pole, it is necessary to find logarithmic (pi — 4) power counting in both the A^, 
i = 1,2, individually. When any virtual eikonal line carries both loop momenta, we can neglect the dependence 
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+ c.c. 



+ c.c. 



FIG. 6: Diagrams that contribute to M2, Eq. (|32p . similar to those studied in Ref. [T(|. These NNLO diagrams describe the 
transition of a pair of momentum P from octet to singlet by the exchange of color with a spectator whose momentum, P3 is 
nearby in phase space. We sum over all cuts of these diagrams that can produce a color singlet quark pair. 



of that line on the softer momentum. We can thus do power counting in the softer loop alone, at fixed values 
of the larger loop momentum. If p so ft < 4 there is at most a single pole. On the other hand, if we find any 
one-loop subdiagram with logarithmic power counting for a given final state, dimensional analysis ensures that 
there is a double pole. 

2. The virtual spectator eikonal lines in the P3 direction (middle lines of the diagrams), contribute A^ 1 , with i 
corresponding to the larger of the loop momenta that flow through the line. 

3. When there is a single active pair eikonal line at the top of a diagram in the amplitude or complex conjugate, it 
gives A~ 2 , with i the larger loop momentum carried by the line. The extra factor of A" 1 relative to the normal 
eikonal line corresponds to a modification of the standard eikonal Feynman rules [l6| due to the integral that 
defines the adjoint eikonal (double) line ending in the field strength vertex, Eq. (|3"2")l . We have encountered this 
effect above in the one-loop calculation, Eq. (pM)) . The three-point crossed vertices, however, each contribute a 
numerator factor of Ai because of the linear dependence on loop momentum in their definition, and four-point 
crossed diagrams are independent of loop momenta. The net effect of the combination of a field strength vertex 
and a single active pair eikonal line is an overall A" 1 when a gluon of momentum fc, attaches to the vertex. 

4. When there are two active pair eikonal lines, as in diagram (V), the diagram represents the sum of two terms, 
in which one of the two eikonal denominators is squared. For example, taking k\ as the gluon connecting the 
field strength vertices, and ki the momentum of the gluon flowing (up) from P3 to the active eikonal in diagram 
(V) , the two active eikonals on the right correspond to the combination 



(P ■ fci) 2 P ■ (fci + fc 2 ) P-h (P ■ (fci + fc 2 )) 2 ' 

We are now ready to discuss the diagrams of Fig. [5] individually. 
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Diagram (IV): Let us begin the application of the above rules with diagram (IV) of Fig. [6l and show that it can 
generate at most a single infrared pole. We must test all independent choices of soft gluon loops. To derive a double 
pole, we must find one loop for which scaling according to Eq. ([33]) at fixed values of the other loop's momentum 
results in a power A^ 4 . 

Again, let k\ be the momentum of the gluon that connect the two crossed vertices, and suppose k\ is the larger 
loop momentum. It must then flow out from the short-distance vertex on the left, and may flow to the short-distance 
vertex on the right either along the double line at the top, or down to the P3 line along the vertical gluon of the 
figure, and from there to the right-hand hard vertex. 

From our observations above, we need only do power counting for the remaining loop, k%, to check for double poles. 
If the vertical gluon is virtual, the &2 loop flows through the right-hand vertex and the vertical gluon line. Power 
counting for this loop depends on how the k\ loop momentum flows. If fcj flows along the top double line, we find 
P2 = 3, while if k\ flows down to P3 along the virtual gluon, P2 = 2. In both cases P2 < 4, and the ki loop cannot 
produce an independent logarithmic divergence. 

The power counting is a bit different if k% describes a gluon in the final state. Since hi is the softer gluon, k\ cannot 
flow (backwards) through the &2 gluon, but must flow only along the double line on top. The power counting for the 
softer line then again gives pi — 2. 

Finally, if ki defines the harder of the two gluon momenta, we readily verify that the k± loop is infrared finite. In 
summary, Diagram (IV) cannot produce a double pole in dimensional regularization. 

In fact, diagram (IV) is the only diagram in Fig. [5] whose individual cuts lack a double pole. We now show, however, 
that for all of the remaining diagrams except for (I), the real parts of the residues of such double poles cancel in the 
sum over allowed cuts (final states) for each diagram. Imaginary double poles, of course, cancel when diagrams are 
combined with their conjugates. We will find a real double infrared pole from diagram (I). 

Diagram (V): Let k\ be the momentum of the gluon connecting the field strength vertices, and ki the momentum 
of the gluon exchanged between the active pair and spectator P3. From the power counting rules described above, 
it is clear that a double pole is produced by taking k\ as the softer loop momentum. There are two final states to 
consider, one for which the gluon ki is virtual, and the other in which it is "real", that is, it appears in the final state. 

Suppose that gluon ki is virtual. Taking the Ai — * limit, the ki integral becomes the same NLO expression as in 
Eq. (p~2|) above (for diagram (V) it is actually the complex conjugate). As in that case, there are two contributions, 
associated with poles in the complex k^ plane at fcf = and P ■ &2 = 0. The former is purely real, the latter purely 
imaginary. When combined with the k\ integral, they both produce l/e 2 contributions in dimensional regularization. 
The imaginary pole, of course, cancels against the complex conjugate diagram. This leaves the real double pole, 
associated with an on-shell gluon &2. This singularit y c ancels against the final state in which the gluon ki is real, 
point-by-point in the remaining momentum integrals [161 ] . 

Diagram (VI): In diagram (VI) double poles arise from a momentum configuration that is complementary to that 
of diagram (V). By analogy to the latter, we choose k^ as the gluon attaching the active pair with P3, and k\ the 
remaining gluon, emitted from the field strength vertex on the left of the diagram. The double pole then arises by 
taking ki as the softer loop momentum. At fixed values of k\, however, the infrared poles arising from the &2 integral 
cancel just as in the NLO case. Again, no double pole survives when diagram (VI) is added to its complex conjugate. 

Diagram (III): In diagram (III), there are actually two relevant loop assignments for which the softer loop has 
logarithmic power counting. To describe them, we denote the momentum flowing out of the field strength vertex on 
the left as k\, and the momentum flowing into the field strength vertex on the right as k^. The momentum of the 
gluon attached to the P3 (middle) line is then &2 — &i (flowing up). The double logarithmic scalings are then Ai ^> A2 
and Ai ~ A2 3> A2 — Ai. In the latter case, the larger momentum flows out of the active pair eikonal and back, and 
in the former case, it flows down to the P3 eikonal line. Both cases correspond to double logarithmic power counting, 
and in the calculation of Ref. [l(| , both regions contributed to a real single pole. 

The calculation of (III) was described in great detail in Ref. [l6| , and the leading, double poles were found to follow 
the same pattern as for diagrams (V) and (VI) above. That is, real double poles cancel between final states in which 
the gluon that attaches to the P3 line is either real or virtual, leaving purely imaginary double poles only. The loop 
integrations that led to this result are the same as in our case, except that the spectator was taken massless, P| = 0. 
Relaxing this condition, it is straightforward to verify that the calculations follow exactly the same pattern as in [l6[ , 
and that diagram (III) does not contribute a double pole after a sum over final states. It is worth noting that most 
of the complications of 16] result from going to the level of the subleading, real single pole. 

Diagram (II): For diagram (II), the pattern is the same as for (V) and (VI): the real double poles are associated 
with on-shell gluons, which cancel, leaving only an imaginary double pole proportional to the result found at NLO in 
Eq. ([30]) . The pattern is again exactly as in the case of a massless eikonal spectator [lj|. 

Diagram (I): We are left with diagram (I) as the only remaining potential source of a real double pole. Diagram 
(I) is associated with final states that include one, two or no gluons. Exactly as in the massless case [16] . the final 
states with one or two gluons cancel against the corresponding contributions from NLO virual diagrams with one or 
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two gluon lines on-shell. (Notice that two-loop virtual corrections do not contribute in the octet-to-singlet transition 
amplitude, M.2-) This real/virtual cancellation, however, leaves the absolute square of the imaginary single-pole term 
of Eq. (|23|) . identified in Fig. [5] Again, this contribution is proportional to P3 , and was hence absent in the massless 
case studied in [l6|. The complete infrared sensitive result thus takes the form 



which exhibits a strong dependence on relative veloctiy. 

We should, of course, emphasize that we have not evaluated the single poles in these diagrams. For massless 
spectators, these poles can be factorized as in Ref. [HI, but in this case we may safely assume that like the double 
poles they will depend on the kinematics of heavy particles in the final state. In any case, having found uncancelled 
double poles, we have already demonstrated the infrared sensitivity of color transfer, independent of the structure of 
the single poles. In the following section, we discuss possible physical implications of these uncancelled poles and the 
consequent infrared sensitivity to color transfer in associated production. 



V. COLOR TRANSFER IN HEAVY QUARKONIUM PRODUCTION 

We have shown that the imaginary 1/e NLO corrections in Eqs. (|23f and (|30|) result in real contributions to the 
octet-to-singlet transition amplitude, M2, at NNLO. As illustrated in Eq. (f5S"|) these corrections are conveniently 
expressed in terms of (3s, Eq. (fT9|) . which measures the relative velocity of the spectator and the pair, as in Eq. (|20|) . 
We are naturally led to suggest that these infrared sensitive corrections tend to increase the cross section for bound 
states, since we have found a new source of pairs with singlet color at small (3s ■ This does not yet give us an estimate 
for its effect on the production cross sections for quarkonia. To make such an estimate, we rely on an analogy to the 
physical picture at the basis of the color octet mechanism as it appears in NRQCD. 



A. Estimating color transfer 



Consider first the standard color octet mechanism. In this case, a (real) infrared pole is associated with the transition 
from a color octet to color singlet pair at NLO through an electric dipole transition [l|. Each such transition is 
associated with a factor of the active pair relative momentum, v in the amplitude, and v 2 in the cross section. We 
can summarize this perturbative result as fTsl [l9j 

d^JcI! ~ da e+e -_ Q Q [Si] (p H ) - ^-v 2 , (36) 

with ph the momentum of the heavy quarkonium, which is identified with the momentum of the active pair. In 
NRQCD, such infrared poles are matched with (equivalently, factorized into) color octet matrix elements, (£%), for 
the final-state quarkonium, H. In this notation, C denotes the orbital angular momentum of the relevant operator. 
For J/tp and similar quarkonia, L = operators give the largest contributions, which requires at least two electric 
dipole transitions. This produces a greater suppression of v A in the cross section, and hence the corresponding matrix 
element, ( 3 S%). We thus make the replacement 

dv e + e -^QQls a ](PH) j ^v 2 -> da e+e -^ QQ[S8] (p H ) ( 3 Sf ) , (37) 

in which the infrared-sensitive correction at NLO is matched to the nonperturbative S'-wave matrix element, even 
though the latter has different scaling in v. 

Now consider the NNLO color transfer cross section in associated production, including the square of the full 
^-dependence at lowest order in v, given in Eq. (|20f . Corresponding to Eq. (|36|) . we have 

pt 1 (1 — (3 2 ) 3 
driller ~ da e+e -^ QQ[Sg]+QI(Ps) (p H ) — a 2 s v 2 .. (38) 

6 Ps 

The kinematic enhancement in this cross section associated with small (3s is very strong whenever (3 < y/v, corre- 
sponding to low relative spectator- active velocities. The quadratic independence here is just the square of the linear 
v in Eq. which is the same as that of an electric dipole transition. As for the color octet mechanism, however, 
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we anticipate that two transitions proportional to the dipole moment (and hence to v) will be necessary to produce a 
color singlet, S-wave quarkonium. In the associated production cross sections at hand, the anticipated v 4 suppression 
that is normally absorbed into the matrix element is compensated, at least in part, by the explicit factor l//?f. In 
effect, for color transfer the expansion in v alone is replaced by an expansion in v//3. 

With our analogy in mind, it seems natural to replace the factor a 2 v 2 /e 2 in (f38|) by the same octet matrix element 
as in (|37| . In making such a replacement, we are assuming that for color transfer it costs the same overall factor of 
ij 4 as in the color octet mechanism to produce an S-wave, color singlet state for the active pair, and also that once 
the active pair is in a color singlet state, it evolves independently of the spectator. We thus summarize our estimate 
of the full cross section for color transfer by 

dal¥ e -_> H+x {p H ) ~ da e+e -^ Q Q [Sl]+Q , {f3s) {p H ) ( 3 S?) 

+ da e+e -^ QQ[Ss]+Q , {Ps) { PH ) &l (1 " Psf ■ (39) 

We emphasize that our arguments give this expression at best a heuristic justification, although we consider it a 
conservative estimate of the color transfer contribution to associated production. We do not rule out the possibility 
that yet higher orders might lead to greater enhancement in /%, but we leave this to further investigation. Another 
approach to estimate color transfer in the velocity-ordered region is based on the observation that the active pair is 
not truly on-shell, but should be treated according to pNRQCD methods [H, 0| • The transverse momentum integrals 
that result in 1/e poles in Eq. (fT3i|) . for example, would then be replaced by logarithms of the bound-state relative 
velocity. In any case, to get a better idea of the role of color transfer, we should compare the sizes of the individual 
singlet and octet hard-scattering functions in Eq. (|39j) , as they appear in the three-particle phase space of the heavy 
quarkonium and the associated pair. We turn now to this disucssion. 



B. Comparison of octet and singlet hard scattering 

The color transfer mechanism is most relevant to the kinematic region where the active heavy quark pair is in close 
vicinity with another heavy (anti)quark. As observed above, color transfer could significantly affect the rate of heavy 
quarkonium associated production, such as e + e~ — > J/ip + cc + X, when the collision energy y/s is not too much 
higher than the mass threshold, 4m, and the cross section is dominated by the region of phase space where one of the 
spectator heavy quarks is close to the active heavy quark pair. In addition, the mechanism should also have a strong 
effect on the size of the input distribution of heavy quark fragmentation functions, e.g., D c ^j/^(z, jio, m), when the 
fragmentation scale fi is close to 3m. In either case, the radiation of light quanta will tend to decrease the effective 
s, and make it easier to produce spectators at low /3. 

At collision energy yfs = 10.6 GeV, as seen by the BELLE and BABAR collaborations, both heavy spectators in 
J/ip associated production can be highly rclativistic if the two pairs carry all the energy. We will find evidence in 
this subsection that even in this case the influence of the color transfer mechanism on J/ip production in association 
with cc could still be significant. 

From Eq. the lowest order inclusive prompt production of a heavy quarkonium with an additional QQ in 
NRQCD is given by 

a ^+e~^H+QQ( S ) = ^e+e--,QQ[S 1 ]+Q'Q'( s ) ( 3S l) 

+ °e + e-^QQ[Ss]+QiQ>( S ) W ) ■ (40) 

Although the direct octet contribution from 3 S to the J/ip production rate at y/s — 10.6 GeV is only about three 
percent of the singlet contribution (20j . the production rate, or coefficient function, for the active heavy quark pair 
is actually much higher for the octet mode than for the singlet, <J e + e -—, C cis 8 ]+c'c' ( s ) ^ ^e+e-~*cc[Si]+c'c' ( s )- The 
octet mode is suppressed because the octet NRQCD matrix element { 3 S%^} is much smaller than the singlet matrix 

element ( 3 <S/^), by about two orders of magnitude @. In the estimate we have given above, Eq. (f39|) . however, the 
contribution via the color transfer mechanism could be very significant if the enhancement factor in [3s compensates 
for the suppression of the octet matrix element relative to singlet. 

The inclusive rate of charmonium production associated with an additional pair of heavy quarks in e + e~ annihilation 
has been studied extensively in the NRQCD formalism at both leading order (LO) 0, Hl| and next-to-leading order 
(NLO) The LO analytic expression for the singlet production of various charmonium states in terms of the 

variable z = 2E c5 /\/~s is available [13, HH, while numerical results for the inclusive rate are available for the NLO 
contribution [l4| • To better understand the contributions to the inclusive rate of prompt J/ip production from different 
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FIG. 7: The integrand, in pb/GeV 3 , of the S3 S4 integration for the production rate of an active heavy quark pair, defined in 
Eq. (142[) . with the pair in a color singlet (left) or a color octet (right) state. 



parts of phase space, we express the rate in terms of a phase space integration over the invariant masses of the active 
heavy quark pair and the heavy spectator (anti)quark: 



s 3 ={P + P 3 ) 2 , Si = (P + P 4 ) 2 



(41) 



with P 2 = (2m) 2 and P 2 = P 2 — m 2 . We denote the total squared collision energy by s = (P + P 3 + Pa) 2 . The 
LO perturbative coefficient functions (or at this order, the partonic cross sections to produce the state QQ[S n ]) in 
Eq. Ignp is given by 



d<J e + e -^QQ[S n ]+Q'Q 



1 Yl x \ Tr l A e+e-^Q*(P/2)Qi(P/2)+Q'Q''P^} e A ( P ) ( 3 * 3 j M 2 



d 3 p 



d 3 p. 



d 3 p 4 



a e 



(2i:f2E QQ {2nf2E QI (2i:f2E Q , 

2 a s{sf 



(2n) i 6 4 (P e+ +P e - -P-P 3 



12s 2 



(-g^W^ (s, s 3 , 54, m)) 9 ($(s, s 3 , s 4 , m)) ds 3 ds 4 , 



(42) 



where n = 1,8, c„ = 2N Cl N 2 - 1, and (3i3j\n) = 8ij/\/W c , V2(T Q ).,; for the singlet and octet contributions, respec- 
tively. The overall factor 1/3 = 1/(2 J + 1) reflects J = 1, the spin of the active QQ pair. In the trace, the projection 
operator for a spin-1 QQ pair is — (7 • P/2 — 771)7** (7 • P/2 + m)/V8m 3 22]. In the second equation in Eq. (I42|) . 
(To = 47ra 2 TO /3s, and e c = 2/3 is the charm quark fractional charge. Finally, explicit forms for the hadronic tensor 
with n = 1, 8 and the phase space constraint 4>(s, S3, S4, m) are given in the Appendix. 

In Fig. [3 we show the integrand of the S3 S4 integration for the production rate a e + e - ^qq^ + qiqi (s) defined in 
Eq. AS)) with an active heavy quark pair in a singlet (left) and octet (right) color state. To generate the figures in Fig. [71 
we used y/s — 10.6 GeV, m = 1.5 GeV. The units for the vertial axes are pb/GeV 3 . After integrating over the phase 

space of S3 and S4, multiplying the same sing let NRQCD matrix element to J/ip, ( 3 5i) = = 1.16 GeV 3 , 

and using the same value of a s used in Ref. [2(j, we obtain the same 148 fb cross section from the direct singlet 
contribution. We also find that the direct octet contribution is about three percent of the singlet contribution. Since 
the octet NRQCD matrix element to J/ip, ( 3 S S ) = (O^d)) = 1.06 x 10~ 2 GeV 3 , is about 100 times smaller than 
corresponding singlet matrix element, it is clear that before multiplying by the NRQCD matrix elements, as shown in 
Fig. [Jj the production rate for a color octet active heavy quark pair is larger than the rate of producing a singlet pair 
everywhere in the phase space including the region when the invariant mass S3 or S4 is small. The octet production 
rate is much higher than the singlet rate when both S3 and 54 are large because the active quark and the antiquark 
of the singlet contribution at this order cannot come from a pair of heavy quarks that originate from either the 
same virtual photon or gluon. With an additional enhancement factor from the color transfer mechanism, the octet 
contribution will also peak in the region where the invariant, S3 or 54, is small. At the same time, the enhancement 
in the color octet coefficient functions away from these regions seems to imply that the cross section for color transfer 
will decay less rapidly than might be expected on the basis of the explicit factors of fts alone in Eq. 
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VI. HIGH ENERGY BEHAVIOR AND FRAGMENTATION 



Equation ([39]) estimates the effect of color transfer on the production cross section as an energy-independent factor 
times a perturbative cross section. In this section, we show why this is a natural assumption, even though the color 
transfer process acts in a restricted region of phase space, specified by Eq. (|2Tj) . It is necessary to check that color 
transfer survives in the high energy limit, and that it can play a role in the fragmentation of heavy flavor at high 
transverse momentum. 



A. Color transfer in fragmentation 

Figure [8] shows a typical fragmentation process for associated production at lowest order, in cut diagram notation. 
To be definite, we discuss fragmentation in leptonic annihilation, as in Fig. [T] The underlying factorization that 
justifies fragmentation analysis extends as well to hadronic scattering, as shown in Ref. [l7l |. 

The shaded circles in Fig. [8] represent as above short-distance functions, now with all propagators off-shell by Q 2 , 
where Q is the center of mass energy. We assume all P 2 = m 2 , taking both heavy flavor masses equal. All propagators 
enclosed within the dashed lines are then off-shell by at least 4m 2 , and remain short-distance contributions from the 
point of view of NRQCD, as in Fig. O 

From our previous considerations, we know that the color transfer process arises from truly infrared gluons, and 
will thus not affect the propagators of lines that are off-shell by the scales m or Q. Rather, the infrared enhancement 
will be associated with factors that depend sensitively on the relative velocities of the external lines. We will now 
show that color transfer is leading power in Q, although it will not appear in logarithmic enhancements associated 
with large relative velocities between P 3 and Pi or P2. 

Let us denote by d,T n (Q^) differential n-particle phase space at total four-momentum Q M . The inclusive associated- 
pair production cross section at total momentum Q 11 = (Q,0) can then be represented as 

da({Pi}) _ 1 f d 3 P 4 



dMQn = ^ I ; dr 3 (Q" - Pt) I M 4 (Q, Pi . . . P 3 ) 



8(2tt) 3 Q 3 



dn 4 J d Sl23 ^ ( Q2 + ? ^q 5123 ) - ™ 2 / dr 3 (Q" - Pi 1 ) I M 4 (Q, P . . . P 3 ) | 2 . 

(43) 



Here and below, we introduce the notation Sij... = (Pj + Pj + . . . ) 2 , while M n is the amplitude for the production of 
n heavy particles. A sum over final-state spins is assumed. The vector P4 is on-shell, with the direction of its spatial 
momentum specified by the angular integral, J dfl 4 . 

At high energy and fixed invariant mass y / si23 of the three-particle (1,2,3) system, the amplitude factorizes into 
a hard function for single-pair production multiplied by an integral that describes the fragmentation of a parent 
(anti)quark into an (anti)quark plus a pair, as in Fig. [5] In the high-energy limit for fixed S123, the corresponding 
squared amplitude integrated over three-particle phase space then factorizes into the squared amplitude for two-particle 
production, times logarithmic integrals characteristic of fragmentation. Suppressing dependence on the relative angles 
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of the Pi, i = 1, 2, 3, the remaining integrals in three-particle phase space appear in the factorized expression as 



/ dr 3 (Q" - Pi 1 ) | M 4 (Q, P ... P 3 ) | 2 ~ | M 2 (P i , Q - Pi) | 2 / ^(Q' 1 - Pf) — 



M 2 {Pi,Q-Pi) 



'123 J 4m 



rf.S 



ds 1 



2:( 



S12 



#(si23 - Sl2 - S23 ~ 3m 2 ) ; 



(44) 



where Q — P4 denotes an on-shell four- vector recoiling against P4. The kinematic dependence in the first line of 
(|4"4")l comes from the denominators of the squared propagators shown in Fig. [51 which give l/sf 2 s? 2 3j multipied by 
momentum factors from the spinor algebra, which give an additional factor in the numerator of Si 2 Si 2 3. A remaining, 
overall factor of 7 • (Q — P4) has been absorbed into the squared two-particle matrix element |M 2 | 2 in the expression. 
In these terms, three-particle phase space reduces to an infrared finite, but logarithmic, integral over the squared 
invariant mass of the (1,2) pair, times a free integral over the squared invariant mass of the (2,3) pair. The theta 
function restriction on S12 and S23 reflects an identity for three-particle phase space: S123 = S12 + S23 + S13 + 3m 2 . 

At lowest order (Fig. [HI for example), the S12 integral in Eq. (|4"4")1 behaves as S123 ln(si 2 3/4m 2 ). This leads in turn to 
an additional logarithmic integral over S123 in Eq. (|43|) . These are standard logarithmic enhancements associated with 
the evolution of a fragmentation function. If we demand, however, that one of the two pairs of Fig. [5J form a heavy 
quarkonium, we will restrict either S23, for the lowest-order color singlet contribution, or s 12 , for the lowest-order 
color octet, to remain close to 4m 2 . The available phase space for the formation of the heavy quarkonium is set by 
the requirement that the active pair should have relative momentum of the order mv. This limitation on phase space 
is taken into account in NRQCD matrix elements, which always include a scaling of order v 3 . In lowest-order singlet 
production, this limitation would be placed on s 2 3, which eliminates both of the logarithmic enhancements identified 
above. For the lowest-order octet mechanism, the restriction can be placed on si 2 , and we may still generate a single 
logarithm from the integral over si 2 3. 

We can now consider infrared-sensitive NNLO color transfer corrections to Fig. [51 In this case, we generate additional 
dependence on the relative velocities of the heavy particles, and hence on the invariants S23 and si 2 . As we have 
seen, this dependence arises from the square of the NLO correction given in Eq. (|23p . and it decreases rapidly with 
increasing relative velocities between the active pair and the spectator. Thus, such a correction suppresses the cross 
section when si 2 3 grows to a multiple of 9m 2 . 

To be specific, recalling the definition of the variable (3s in Eqs. (jTUJ) and (fT9"]) . and neglecting v compared to 0g, 
we find an additional si 2 3-dependence in the integrand of Eq. (J44J) , 



(i - Bl? 



a dsi2 at 

ds 2 3 0(S12 3 

Sl2 



S12 - S23 

64m 6 



3m 2 ) 



is 



123 



9m 2 )(si23 — 5m 2 )(si23 — m 2 ) 



ds 



ds 



12 
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S12 



(S123 - su - s 2 3 - 3m 2 ) 



(45) 



where the lower limit of the si 2 3 integral will be set by (3s > v, and where there is a rapid decrease for S123 S> m 2 . As 
with the color octet mechanism for this lowest-order process, the (1,2) pair forms the bound state, and we anticipate 
an overall suppression in the S12 integral of v 3 on dimensional grounds. Although non-logarithmic, the S123 integral 
remains leading power in Q, with no further suppression by powers of v. 



B. Light quarks and polarization 



By referring to Fig.[5J we can make a very simple observation relevant to the comparison of the associated production 
cross section for heavy quarkonia with heavy flavors compared to light flavors only. The leading power behavior for 
the production of single heavy quark pairs in close proximity in phase space is due entirely to fragmentation. The 
color transfer contribution for fragmentation diagrams as shown in the figure is present only for associated production, 
since light quarks will unavoidably have much larger relative velocities, Ps- At lowest order, then, we would expect 
the fragmentation of a light quark to a heavy quarkonium to proceed only through the color octet matrix element. 
For leptonic annihilation in particular, where gluon fragmentation is yet higher order in a s , the color transfer process 
affects heavy-flavor associated production only. 

Related considerations .apply to the prediction of transverse polarizations for vector heavy quarkonia produced from 
gluon fragmentation [Tol. [T1LT23II . It is easy to check that the color transfer process does not give such a prediction. The 
color transfer itself respects the spin of the pair, which reflects the polarization of a collinear gluon in the fragmentation 
process (the gluon of invariant mass S12 in Fig. [5J. Indeed, at high energies, the polarization of this gluon naturally 
includes large longitudinal components. A further study of this issue is clearly in order. 
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VII. SUMMARY AND CONCLUSION 

The color transfer process that we have described above provides a new viewpoint on the hadronization of heavy 
quarks in associated production. In color transfer, the pair produced from a single gluon can be transformed from 
octet to singlet representation by the field of an open (anti)quark that is at sufficiently low relative velocity. We can 
picture it as a process that catalyzes quarkonium production when another heavy quark is nearby in phase space. 
Experimental confirmation of such a process could provide a significant tool to study dynamical processes in quantum 
chromodynamics. Color transfer complements the standard picture of hadronization through string breaking, in which 
color neutral pairs are formed from quarks and antiquarks that originate from separate virtual gluons [25| . The latter 
is leading power in the number of colors, and color transfer is hence nonleading in this expansion. 

We have shown that in perturbation theory, color transfer appears first at NNLO. At this level, it is positive, but 
not infrared safe. Although nonperturbative it cannot be described in terms of the matrix elements of NRQCD in 
general. Without experimental input, we are not yet in the position to give a realistic estimate of its importance. 

Color transfer could be part of the explanation of well-known anomalies in comparisons of data to NRQCD pre- 
dictions [Tol . [l3j . and it may be possible to confirm its role through a number of qualitative predictions that flow 
naturally from the perturbative reasoning above. 

First, we would expect color transfer to produce phase space distributions that are peaked at small relative invari- 
ant mass between the heavy quarkonium and closest open heavy flavor in associated production. For example, an 
experimental signal of this effect would be an enhancement of J/ip production relative to NRQCD estimates at low 
values of the invariant masses of J/tp — D and J/tp — D pairs. At high energies, a similar effect should be found both 
in T — D and T — B systems. Second, if color transfer enhances associated production at large transverse momenta, 
it could also help explain the observed polarizations of heavy quarkonia. 

A Dalitz plot analysis of quarkonium/open flavor final states at B factory energies could show enhancements on 
the low pair mass corners. At the same time, our analysis suggests that in the interior of such a Dalitz plot, an 
NRQCD analysis based on the color singlet mechanism should succeed. In general, semi-inclusive measurements of 
the distribution of heavy flavor in final states could even enable us to probe the time evolution of flavor, momentum 
and spin in the formation of the final state. This would complement the inclusive viewpoint which is built into the 
important NRQCD calculations of Refs. [HI, 03] ■ 

Color transfer is distinct from, yet in some ways analogous to, the usual NRQCD color octet mechanism, which relies 
on soft gluon radiation. It is kinematically dependent on soft gluon exchange with other particles in the final state, 
and hence is not "universal" in the sense of NRQCD. On the other hand, the effect decays with relative transverse 
momentum, and hence does not spoil high-pT factorization, as described, for example, in [16j |. Similarly, it does not 
occur at all through exchanges with massless particles, and hence is specific to associated production with heavy 
flavor. 

In conclusion, although many issues remain to be studied, the color transfer mechanism suggests a number of 
phenomenological signals, which should make it possible to test its relevance to associated production. If it does pass 
these experimental tests, it may offer insight into the dynamics of color in hadronization. 
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APPENDIX A: LOWEST ORDER MATRIX ELEMENTS FOR HEAVY QUARKONIUM ASSOCIATE 

PRODUCTION 

In this appendix, we provide the leading order squared matrix elements, —g^W^ that appear in Eq. (|42[) . for 
producing an active pair of heavy quark and antiquark of mass to, in association with another heavy quark pair of 
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the same mass. For the active pair in a singlet color state, n = 1, we have 



„ = <il [ 2 

27m \ ( m 2 _ S3 ) 2 ( m 2 _ S4 )2 (_ 2m 2 + S3 + fi4 _ 2s )2 



132m i0 + 2(136s - 101(s 3 + s 4 ))m s 



+4 (7s 3 + 108s 4 s 3 + 7s 2 4 + 240s 2 - 58(s 3 + s 4 )s) m 6 

- (31s| + 149s 4 s 3 + 149s 4 s 3 + 31s 3 4 - 380s 3 + 316(s 3 + s 4 )s 2 - 9 (9s 2 . - 2s 4 s 3 + 9s^) s) m 4 
+ (s 4 + 56s 4 s| + 98s^s^ + 56s 3 4 s 3 + s\ + 32s 4 - 60(s 3 + s 4 )s 3 

+12 (3sj - 2s 4 s 3 + 3sf) s 2 - (s 3 + s 4 ) (s§ - 10s 4 s 3 + sf) s) m 2 

+s 3 s 4 (-9S3 - 25s 4 s 2 - 25s 2 s 3 - 9s 4 - 4s 3 + 4(s 3 + s 4 )s 2 + (9s 2 . + 14s 4 s 3 + 9s 2 ) s) 
f 64m 6 + (41s 3 + 9s 4 + 47s)m 4 + 2s 3 (s - 5(s 3 + s 4 ))m 2 + s 2 (s 3 + s 4 - s) 



(m 2 - s 3 ) 4 

2 (m 2 - s 3 ) r R 1 

+ 5^ ^ 257m 6 + 3 2s 3 - 7s 4 + 54s)m 4 

—2m 2 + s 3 + s 4 — 2s 

- (7sg + 36s 4 s 3 + 10s 2 - 16s 2 + 2(s 3 + 7s 4 )s) m 2 + s 3 s 4 (s 3 + 2s 4 - 2s)] 

4 ( 771" - S'i) ^ 

+ i ^ 2 [!60m 6 + 3(5s 3 - 15s 4 + 37s)m 4 

(—2m 2 + S3 + s 4 — 2s) 

-2 (5s 2 + 14s 4 s 3 - 2ss 3 + ls\ - 4s 2 + 5s 4 s) m 2 + (2s 2 + 4s 4 s 3 + 3s 2 ) (s 3 + s 4 - s)] (1) 
64m 6 + (9s 3 + 41s 4 + 47s)m 4 + 2s 4 (s - 5(s 3 + s 4 ))m 2 + sl(s 3 + s 4 - s) 



4 

m 



(m 2 — s 4 ) 4 

+ -2m 2 t 2 3 + S s 4 ) -2s [ 25?m6 + ^ + 27a > - 21fl »> 

- (10s 2 + 36s 4 s 3 + 7s 4 - 16s 2 + 2(7s 3 + s 4 )s) m 2 + s 3 s 4 (2s 3 + s 4 - 2s)] 

+ '™ 2 [!60m 6 + (-45s 3 + 15s 4 + llls)m 4 

(—2m 2 + S3 + s 4 — 2s) 

-2 (7s 2 + 14s 4 s 3 + 5ss 3 + 5s 2 - 4s 2 - 2s 4 s) m 2 + (3s 2 + 4s 4 s 3 + 2s 2 ) (s 3 + s 4 - s)] 
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For the pair in an octet color state, n = 8, we have 

1 



9m (to 2 — s 3 ) (m 2 — S4) 



^ 1594m 14 - 4588s 3 m 12 - 4588s 4 to 12 + 4285s 2 to 10 + 4285s; 
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+13670s 3 s 4 m 10 - 1961s^m 8 - 1961s^m 8 - 12723s 3 s 2 m 8 - 12723s 3 i s 4 to 8 + 553s 4 m 6 + 553s 4 to 6 
+5699s 3 s|to 6 + 10316s 2 s^m 6 + 5699s^s 4 m 6 - 75s 3 'to 4 - 75s|to 4 - 1417s 3 s 4 to 4 - 3760s 3 ! s 4 ; to 4 
-3760s|s 4 to 4 - 1417s 3 s 4 to 4 + 127s 3 s|m 2 + 791s 3 s 4 m 2 + 846s^s|m 2 + 791s 3 s 4 to 2 + 127s 5 3 s 4 m 2 
-61sjsl - 137s^s^ - 137s^s| - 61s|s^ 



4 (m 2 — s 3 ) (to 2 — s 4 ) ■ 



-2to 2 + s 3 + s 4 — 2s 



478m 12 - 1130(s 3 + s 4 )to 10 + (619s 2 . 



3084s 4 s 3 + 619s 2 ) to 8 



-4(s 3 + s 4 ) (74s 2 . + 383s 4 s 3 + 74s|) m 6 + 2 (84s^ + 256s 4 s^ + 441s 2 s 2 ! + 256sls 3 + 84s 4 ) m 4 
-(s 3 + s 4 ) (l9s| + 110s 4 sJ + iOsjsj + 110s|s 3 + 19s 4 ) m 2 

+s 3 s 4 (7s 3 + 21s 4 s 3 + 6s 2 s 3 + 21s 4 s 3 + 7s|) 
2 o \ 2 /™2 



(2) 



4 (to' — S 3 j (TO' — s 4 



-2to 2 



s 3 + s 4 — 2s) 



244to 10 - 458(s 3 + s 4 )to 8 + 4 (17s 2 , + 304s 4 s 3 + 17sf) to 6 



-2(s 3 + s 4 ) (57s 2 + 136s 4 s 3 4- 57s 2 ) to 4 + 2 (l9s 3 
-(s 3 + s 4 ) (3s| + 8S4S3 + \2s\sl + 8s|s 3 + 34 
2 (7to 2 - 3s 3 ) (m 2 



52s 4 s 3 3 + 52s 2 4 s 2 3 + 52sis 3 + 19s 4 ) to 2 



S3 



) (7to 2 - 3s 4 ) (to 2 - s 4 ) s 2 9s|s| (6to 4 + 2(s 3 + s 4 )to 2 + s 2 + s|) 



2s 



252to 14 - 677(s 3 + s 4 )r 



(SOlsj + 2438s 4 s 3 + 301s^) to 10 



+ 0s + s 4 ) (173s 2 . - 2185s 4 s 3 + 173s 



TO- 



-97s^ + 410s 4 s^ + 1946s^s^ + 410s^s 3 - 97s^) to' 



+s 3 s 4 (s 3 + s 4 ) (33s 2 - 581s 4 s 3 + 33s 2 ) to 4 + s 2 s 2 (25s 3 + 166s 4 s 3 + 25s 2 ) to 2 - 9s 3 s|(s 3 + s 4 ) 

The invariants in these expressions are defined as s = (P + P3 + P4) 2 , s 3 = (P + P 3 ) 2 , and s 4 = (P + P4) 2 , respectively. 
The phase space integration of ds 3 ds 4 in Eq. (142|) is given by the condition, $(s, s 3 , s 4 , to) > 0, with the function 



<I>(s, s 3 , s 4 , to) = s 3 s 4 (6to 2 — s 3 — s 4 + s) + 3(s — to 2 )(s 3 + s 4 )to 2 — (s + 6to 2 )(4s + to 2 )to 2 + 2(9s + 4to 2 )to 4 . (3) 
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